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The major challenge in cancer therapy is to efficiently translocate drug molecules into cancer tumors 
without doing any damage to healthy tissues. Since there exist pH gradients between tumor and normal 
tissues, pH- sensitive materials may have great potential to overcome such challenge. Here, we report one 
new type of pH-responsive drug delivery system where pH-sensitive polymers are introduced to control the 
cellular uptake of nanoparticles under different pH environments through dissipative particle dynamics 
simulations. Interestingly, the behavior of cellular uptake of nanoparticles here exhibits "smart" 
pH-responsive properties: for lower and higher pH, the nanoparticles can be taken up by cell membranes, 
while for pH in middle range, the endocytosis is blocked. Further, it is found that receptor-ligand 
interactions as well as surface charge property of nanoparticles and membranes can also have important 
impacts on the endocytosis. The present study may give some significant insights into future 
stimulus-responsive medical materials design. 

Efficient delivery of nanoparticles into cell interiors is of great importance in biomedicine^'^. For example, 
how to efficiently transport drug molecules into cancer tumors without doing any damage to healthy tissues 
is still a huge challenge in cancer therapyl Recently, the pH- responsive materials which can exhibit a sharp 
change in physicochemical properties under different pH environments have shown a great potential to realize 
such functions^"^, because some solid tumor tissues form a mildly acidic environment (pH — 6.5), lower than that 
in the normal tissues and blood (pH ~ 7.4)^"^. 

However, there is still little knowledge about the cellular uptake of nanoparticles, especially for how nano- 
particles translocate across cell membranes under different external pH. Although some previous experimental 
studies have shown that the change of pH may be utilized to transport nanoparticles into targeted tissues^"^°, it is 
rather difficult to systematically probe and visualize the endocytosis process under various conditions, due to 
available experimental technologies. Furthermore, the in vivo and in vitro efficiencies are still low because there 
exist lots of barriers for tumor drug delivery^ \ It is therefore essential to have a thorough understanding of the 
detailed molecular mechanism of the translocation process under different pH environments. Moreover, since 
the physicochemical properties of nanoparticle and the ligands on its surface may have important impacts on the 
cellular uptake of nanomaterials^^"^^, can we design any new types of stimulus materials to make the best use of the 
variable properties of nanomaterials and well control the endocytosis process under different environments? 

In this study, we undergo the first computational study to design such type of pH-sensitive nanomaterials (i.e., 
nanoparticle-polymers complexes) by using dissipative particle dynamics (DPD) simulations Differing from 
the main methods in experiments where nanoparticle (or/and its covalent ligands) may show different charge 
properties under different pH environments^"^, the key idea here is to use external pH-sensitive polymers to well 
control the cellular uptake of nanoparticles. As we will show below, the uptake behaviors here can show multiply 
pH -responsive, which may have some advantages over previous designed materials in real applications. Further, 
the effect of properties of nanoparticles, ligands and membranes on the cellular uptake will also be studied. 

Results 

Idea of pH-sensitive nanoparticle-polymers complex design. Figure 1 shows the coarse-grained models of 
different components in our simulations. The nanoparticle-polymers complex (NPC) is composed of one 
nanoparticle and some pH-sensitive polymers (Fig. la), and can be spontaneously formed due to electrostatic 
interactions when placing the nanoparticle in the polymer solutions (see Supplementary Fig. SI a). The 
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Figure 1 | Schematic illustration of the models in the simulations, (a) Snapshot of the nanoparticle-polymers complex (NPC, i.e., the nanoparticle with 
pH-sensitive polymers assembling onto its surface) ; (b) Snapshot of membranes and architecture of lipids and receptors. Green bead represents lipid head 
containing + e, purple bead is lipid head containing -e, while blue bead stands for lipid head with no charge; the orange bead represents lipid tail, the red 
bead stands for receptor head bead, yellow bead is particle bead, and the polymers are formed of the cyan beads (with charge -e) and pink beads. 



nanoparticle is fabricated by arranging hydrophilic DPD beads (P) on 
a fee lattice with lattice constant a = 0.40 nm into a desired geometry 
shape and volume, and all beads comprising a nanoparticle move as a 
rigid body^^'^^. The surface beads on the particle are treated as ligand 
ones^^'^^. Each surface bead carries a charge of +e (the surface charge 
density is set as ap = 2.5 e/nm^ unless otherwise stated), and its 
ionization degree is kept 100% in the simulations (some chemical 
group like l,8-Bis(dimethylamino)naphthalene^'* may be used to 
reach such demand when pH is not very large). The pH- sensitive 
polymer is composed of twelve hydrophilic beads (C), and some 
beads may carry a charge of — e, which depends on the system's pH 
value (for real applications, many chemical groups like carboxyl group 
can be utilized in the pH-sensitive polymer design^^). For the 
simplicity, we use the relationship between the fraction (0) of 
charged polymer segments and external pH as follows^^'^^: 



1-0 



-\-pKa, where pK^ is the acidity constant of the 



polymer. Note that the charged monomers are randomly generated, 
but with a fixed charge number N per polymer {N = 120 in a fixed 
polymer length of 12) during the simulations^^. By doing some simple 

transformation, we have N = — r-. From this equation, we 

can easily find that N increases with the increase of pH. 

Each amphiphilic lipid consists of a headgroup containing four 
connected hydrophilic beads (H) and two tails with respective three 
hydrophobic beads (T) (Fig. lhy^~^\ The first head bead carries a 
charge of +e, while the second head bead carries a charge of — e; the 
remaining two beads are uncharged^^. Particularly, when modeling 
the negatively charged lipids, non-charged hydrophilic bead is used 
to take place of the first positive charged bead in lipid molecule. The 
receptor molecule has the same conformation of lipid mole- 
cules^^'^^'^^ but its first two head beads (R) are uncharged and can 
interact with the ligand bead (L) via soft Lennard- Jones (LJ) poten- 
tials^^. When lipids and receptors are immersed in the water, they can 
form a stable membrane (see Supplementary Fig. Sib for the density 
profile of all types of beads). Additionally, the percent of the recep- 
tors in the membrane is set to be 50%^ 
radius is fixed as 4 nm. 



and the nanoparticle 



Endoqiiosis of nanoparticle-polymers complex under different 
external pH. As shown in Fig. 2, there exist big differences in the 



endocytosis process under different pH values — the endocytosis of 
nanoparticle will be blocked by the associated polymers on its surface 
under middle-range pH, whereas the nanoparticle can be fully 
engulfed by the membrane under low and high pH. 

When pH is low (pH < pK^ - 0.48, i.e., N < 3), the number of 
ionized monomers is small. As a result, the adsorption of polymers 
on the particle surface is so weak that the polymers cannot prevent 
the strong attraction between the receptors on the membrane and 
ligands on the nanoparticle surface. Therefore, the polymers will 
detach from the nanoparticle surface with the increase of the engulf- 
ment degree. Finally, the particle can be totally engulfed by the mem- 
brane, with no polymers on its surface (Fig. 2a). 

When pH is in middle range (pK^ - 0.30 ^pH^ pK^ + 0.30, i.e., 
4 < AT < 8), there are several monomers ionized in each polymer, 
so the polymers can absorb on the nanoparticle surface stably. 
Additionally, since the electrostatic repulsion between adsorbed 
polymers is not sufficiently strong, the number of adsorbed polymers 
can be very large (see Supplementary Fig. S2). As a result, the nano- 
particle is totally surrounded by the polymers. Given that the 
receptor-ligand interaction is short- ranged, the receptors on the 
membrane cannot interact with the ligands on the particle surface. 
Therefore, in this situation the NPC will just weakly attach to or get 
away from the membrane (Fig. 2b). 

When pH is high {pH > pK^ + 0.48, i.e., N > 9), the number of 
ionized monomers becomes large, which makes the polymers 
strongly absorbed onto the nanoparticle surface. However, because 
of the strong electrostatic repulsion between charged polymer mono- 
mers themselves, the number of adsorbed polymers is small. As a 
result, the polymers cannot totally prevent the nanoparticle from 
receptors so that the nanoparticle can also attach onto the membrane 
and then be engulfed by cell membranes. However, differing from the 
former case of low pH, there still exist some polymers on the nano- 
particle surface after it is totally engulfed (Fig. 2c). This is because in 
this case the electrostatic interaction becomes strong (which is com- 
parable to the receptor-ligand interaction) so that the polymers are 
very hard to detach from the nanoparticle surface. 

To obtain further insight of the three different endocytosis 
processes, we calculate the potential of mean force (PMF) of them 
along the membrane-normal direction by using steered molecular 
dynamics (SMD)^^'^^. Figure 3a shows the PMF as a function of the 
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Figure 2 | Time sequence of the snapshots of interactions between the nanoparticle and membranes in the presence of pH-sensitive polymers under 
different ionized degrees, (a) N = 1 (i.e., pH = pK^ - 1.04); (b) N = 6 (i.e., pH = pKj; (c) N = 11 (i.e., pH = pK^ + 1.04). 



distance (Z) of particle from the central plane of the bilayer when N 
= 1,6, and 11, corresponding to three different engulfment modes 
shown in Figs. 3b-d, respectively. For N = 1, as the nanoparticle 
approaches the membrane, the PMF monotonously decreases and 
there is no energy barrier in the PMF curve in this case, indicating 
that the endocytosis of nanoparticle is spontaneous. While for N = 6, 
the PMF monotonously increases with the decrease of the distance 
(Z) between nanoparticle and membranes, so the engulfment of 
particles is impossible in this case. When N is 11, the PMF firstly 




Figure 3 | Comparison of the three different interaction modes. 

(a) Potential of mean force (PMF) of nanoparticles as a function of its 
distance (Z) from the central plane of the lipid bilayer in three different 
modes, (b)-(d) show the snapshots of final equilibrium of interactions 
between NPC and membranes under different external pH: (b) N = 1, 
(c) N = 6, and (d) N = 11. Error bars are obtained by taking the standard 
derivation of ten independent simulations. In order to distinguish the 
three modes, we name them as total engulfment type 1 (TEl), blocked 
engulfment (BE), and total engulfment type 2 (TE2), respectively. 



increases with the decrease of Z, which is due to small deformation of 
membranes and weak repulsive interactions between polymers and 
membranes. Nevertheless, when the NPC is closer to the membrane, 
the PMF decreases because of the strong attractive receptor-ligand 
interactions. The energy barrier in this case is just about 5 /c^T, and 
could be overcome only by thermal motions, therefore the NPC may 
also be engulfed by membranes. 

Effect of receptor-ligand interaction as well as charge property of 
nanoparticle and lipid membrane on the endocytosis. To deeply 
understand the mechanism of the endocytosis of nanoparticles in the 
presence of pH-sensitive polymers, we further study other factors on 
the nanoparticle internalization. As discussed above, there may exist 
three different phases (i.e., TEl, BE, and TE2) in the interactions 
between nanoparticles and membranes. Additionally, another two 
phases including weak absorption (WA) and partially engulfment 
(PE) can also occur when the receptor-ligand interaction is not 
strong enough according to previous studies^^'^^. As a result, here 
there are five different phases. Note that the difference between 
TEl and TE2 is whether there exist polymers on the nanoparticle 
surface. 

As shown in Fig. 4a, when the receptor-ligand interaction strength 
is weak, there will be no obvious engulfment of the NPC. With the 
increase of the receptor-ligand interaction strength, the NPC will just 
be partially engulfed by membranes and when the strength is strong 
enough, the NPC can be totally engulfed under both low and high pH 
conditions. However, when AT is in the range of 4-8 (i.e., pH is in [pKa 
— 03,pKa + 0.3] ), the engulfment will be blocked because the ligand 
can not interact with receptors, even for stronger receptor-ligand 
interaction. Therefore, the increase of interaction strength in this 
case has little effect on the nanoparticle internalization. 

Besides the interaction strength, the range of the receptor-ligand 
interaction is also a characteristic quantity^^'^^, which may have sig- 
nificant effects on the cellular internalization of nanoparticles. As 
shown in Fig. 4b, with the increase of the interaction range, the 
change of the phases in low and high pH is not obvious except that 
the phase (N = 2, Vcut — 2.0 nm) changes from TEl to TE2 (the 
reason is that when the interaction range is large, the ligand on the 
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Figure 4 | Phase diagrams describing the equihbrium state of the particle-membrane systems, (a) the receptor-Kgand interaction strength verse the 
ionized degree N, where r^wf = 1.0 nm and a p = 2.5 e/nrn^; (b) the receptor-ligand interaction range verse the ionized degree N, where e = lO/CfiTandcrp = 
2.5 elnm^-y (c) the surface charge density verse the ionized degree N, where Ycut = 1-0 nm and e= lO/c^T. In (a)-(c), the membrane is neutrally charged, 
(d) Negative lipid fraction (i.e., the number of negative lipids/the total number of lipids and receptors) verse the ionized degree N, where e= lO/c^T, 
Gp= 2.5 elnrn^, and Ycut =1-0 nm. 



nanoparticle surface can still interact with receptors on the mem- 
brane even if some polymers attach to the nanoparticle). 
Importantly, the NPC can be engulfed by membranes when the 
ionized number is in the range of 4-8 because the polymers on the 
particle fail to protect the "longer" range interactions between 
ligands and receptors. Here the final equilibrium is also TE2, based 
on the same reason discussed above. In this sense, the receptor-ligand 
interaction range may be more important for cellular uptake of 
nanoparticles than the interaction strength. 

Further, we examine the effect of surface charge density of nano- 
particle on the internalization. As shown in Fig. 4c, when the surface 
charge density is zero, the polymers do not attach onto the nanopar- 
ticle due to lack of electrostatic attractive interactions. As a result, the 
engulfment of nanoparticle will not be affected by polymers, and the 
nanoparticle can be totally engulfed with no polymers on its surface 
(TEl). With the increase of surface charge density, the phases in high 
pH case will begin to change from TEl to TE2, because in that case 
the polymers carrying more charges will not detach from the nano- 
particle surface during the engulfment. As the charge density 
increases, more polymers will absorb onto the nanoparticle surface 
so that the engulfment of nanoparticle will be blocked. Even if the 
ionized number AT is small (e.g., N = 2, ap = 3.75 elnrn^), the block 
could also happen when the charge density is large enough. In addi- 
tion, in the case of high pH and high surface charge density, since the 
electrostatic repulsion between highly charged polymers become less 
important (compared to the electrostatic attraction between nano- 
particle and polymers), more polymers can attach onto the nanopar- 
ticle surface to prevent the interaction between nanoparticles and 
membranes. As a result, under that condition (e.g., N = 10, cTp = 
3.75 e/nm^), the engulfment will also be blocked. 

Since there are more anionic molecules (e.g., sialic acid) abundant on 
the surface of cancer cells as compared to normal cells^^, it is very 
important to study the effect of membrane charge property on the 
endocytosis. For the sake of the simplicity, here we just use the nega- 
tively charged lipids to model the anionic molecules in the membrane^^. 



When N < 2, the zeta potential(C)^°'^' of NPC is above 40 mV (see 
Supplementary Fig. S3), so there wiU exist the attractive electrostatic 
interaction between the NPC and negative membranes, which can help 
the cellular uptake process. Thus the final state under low pH will stiU be 
total engulfment. With the increase of N, C decreases largely (when N = 
4, C is about 5 mV, and when N > 6, C is below —30 mV, see 
Supplementary Fig. S3). As a result, the attractive electrostatic inter- 
action between the NPC and membranes becomes very weak when N 
= 4 so that it will not affect the endocytosis process. While when N > 6, 
the electrostatic interaction becomes repulsive, which will greatly pre- 
vent the endocytosis, so the final state in middle pH stiU is blocked 
engulfment. Generally, as shown in Fig. 4d, the negative membrane has 
little impact on the phase behaviors under low and middle pH. On the 
contrary, since the NPC under high pH is negative-charged, when the 
membrane carries more negative charges, the repulsive electrostatic 
interactions will become very strong, which can induce the failure of 
the adsorption of the NPC onto cell membranes. Therefore, the endo- 
cytosis wiU be blocked in this situation (Fig. 4d). 

Discussion 

In the present study, we report one new type of pH-responsive drug 
delivery system with the addition of pH- sensitive polymers for con- 
trolling the cellular uptake of nanoparticles under different pH envir- 
onments. Importantly, as summarized in Figs. 3b-d, the endocytosis 
process shows triple-pH-responsive, i.e., under lower and higher pH 
conditions, the particles can be engulfed by cell membranes, while in 
the middle range of pH, the endocytosis process is blocked. This 
triple-pH-responsive property is quite different from those reported 
by previous studies^"^^ and may have some advantages over them in 
multi-pH-responsive process. For example, the designed NPC with 
triple-pH-sensitive property may have the capability to distinguish 
the multi tissues of various pH values. It is known that there may exist 
pH gradient between different tumor types^^ (e.g., pH in sarcomas is 
about 7.0 while in astrocytomas it is about 6.5). In this case, we can 
choose some type of polymers (e.g., PAMAM-PEG-mannose^^) with 
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pKa of about 7.0. When the environmental pH is about 6.5 (AT ~ 2 - 
3) or 7.5(N ~ 9 - 10), the nanoparticle (e.g., quantum dots) can be 
well translocated into cell interiors, whereas in the case of pH of 
about 7.0 (AT ~ 6), fewer nanoparticles are taken up by cells. 
Therefore, there exist big differences of fluorescence intensity, which 
can distinguish not only the normal tissues from cancer tumors but 
also different tumors. 

Furthermore, the NPC can be also used in drug delivery for cancer 
therapy, where the polymer with pK^ of about 7.5 could be chosen 
(e.g., PLGA-b-PBLC). When pH is about 6.5 (iV ~ 1, typical cancer 
cell environment), the cellular uptake of particles can be very high, 
whereas the uptake will become low under pH of about 7.5 (N ~ 6) 
for normal cell, which can make the best use of drugs. Moreover, 
since here the uptake is low in normal cells, this may also decrease its 
potential toxicity^^. Thus we believe that the high targeted efficiency 
as well as the low nanotoxicity will make it a good type of drug 
carriers for cancer therapy. 

Finally, we point out the feasibility and application of our reported 
drug delivery system in real experiments. To some extent, the present 
drug delivery system may be just suitable for in vitro experiment. 
While for in vivo experiment, since there exist other charged macro - 
molecules^^ which could also affect the cellular uptake, it becomes 
very complicated in that situation and is beyond the present system. 
Nevertheless, if we treat the charged polymers in our system as 
external charged macromolecules (e.g., proteins), the present study 
can give some useful insights into another important question in 
drug delivery — how do external macromolecules affect the cellular 
uptake of nanomaterials? Recently, an experimental study has illu- 
strated that nanoparticles may lose its targeting ability when placed 
in a complex biological environment^^. Here we may give a quant- 
itative theoretical explanation for it, i.e., the adsorption of external 
macromolecules on nanoparticle surface can suppress the specific 
interactions between receptors and ligands. Generally, our results 
show that the cellular uptake of nanoparticles can be well controlled 
by the properties of nanoparticles, pH-sensitive polymers, cell mem- 
branes, and external environment. Therefore, our study can help 
better understand the physical mechanism of receptor- mediated 
endocytosis of pH-sensitive nanomaterials, and may provide signifi- 
cant ideas for engineering new types of stimulus-responsive materi- 
als in biomedicine. 

Methods 

We briefly summarize the simulation method used in this work and the details can be 
found in Supplementary Information. Dissipative particle dynamics (DPD) is a 
coarse-grained simulation technique with hydrodynamic interaction^''. The dynamics 
of the elementary units which are so-called DPD beads, is governed by Newton's 
equation of motion. Typically, there are three types of pairwise forces in DPD, i.e., the 
conservative force, dissipative force, and random force. Here, in order to include 
electrostatic interaction between charged beads, the Coulomb force is incorporated 
into our DPD simulations^"'^^. Further, a modified LJ potential is introduced to mimic 
the receptor-ligand interactions^^. Additionally, we also use a harmonic bond to 
ensure the integrality of lipids and polymers, and a three-body potential to depict the 
rigidity of lipid tails^*^ ''^. All simulations are performed in the NVT ensembles using 
the velocity- Verlet integration algorithm. The size of the simulation box is 65rc X 65rc 
X 40rc with the number density of p = . The integration time step At = 36 ps and 
each simulation time is at least 18.0 fis. 
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